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ABSTRACT 

We present a simultaneous detection of gravitational magnification and dust reddening 
effects due to galactic halos and large-scale structure. The measurement is based on 
correlating the brightness of '^85,000 quasars at z > 1 with the position of 20 million 
galaxies at z ~ 0.3 derived from the Sloan Digital Sky Survey and is used to constrain 
the galaxy-mass and galaxy-dust correlation functions up to cosmological scales. 

The presence of dust is detected from 20 kpc to several Mpc, and we find its 
projected density to follow: Sdust ^ 0^^-^, a distribution similar to mass. The amount 
of dust in galactic halos is found to be comparable to that in disks. On large scales 
its wavelength dependence is described by Ry — 3.9 ± 2.6, consistent with interstellar 
dust. We estimate the resulting opacity of the Universe as a function of redshift and 
find {Ay) ^ 0.03 mag up to z = 0.5. This, in turn, implies a cosmic dust density of 
J7dust — 5 X 10~^, roughly half of which comes from dust in halos of ~ L* galaxies. 

We present magnification measurements, corrected for dust extinction, from which 
the galaxy-mass correlation function is inferred. The mean mass profile around galaxies 
is found to be S ^ 30 (6'/l')^°-^ /iMqpc"^ up to a radius of 10 Mpc, in agreement 
with gravitational shear estimates. 

Key words: dust - extinction, reddening - dark matter - magnification ~ large-scale 
structures - quasars - galaxies 



1 INTRODUCTION 

Light rays from distant sources carry unique informa- 
tion about the matter and gravitational potential along 
the line-of-sight. A well-known example is the signature 
of intervening gas clouds imprinted into spectra of back- 
ground sources via absorption lines. Mass concentrations 
located along the path of photons can also induce gravi- 
tational lensing effects . Background s o urces can be mag- 
nified, as detected by IScranton et al.l |200^, and galaxy 
shapes can be dist orted as m easured through galaxy- 
galaxy lensing (cf. Hoekstra et al. 2002 , Sheldon ct al.. 
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Parker et al.. ,20071 ') and cosmic shear (e.g. 
Measuring these effects has become a powerful tool for prob- 
ing the mass distribution in the Universe. 

In addition, dust extinction effects are expected to occur 
as radiation pressure from stars and galactic winds triggered 
by supernovae are expected to expel some amount of dust 
from galactic disks into the intergalactic medium ijAguirrd 



1 19991 : [ Bianchi fc Ferrarallioosi '). Detecting dust reddening by 
galaxy halos would provide us with useful information on 
the life cycles of dust particles as well as characterize the 
opacity of the Universe. In practice, detecting such an effect 
is made difficult by the requirement to measure brightness 
and/or color changes at a sub-percent level on 100 kpc scales 
around galax ies. A first attempt to find dust in galactic halos 
was made bv lZaritskvl l| 19941 ) who reported a 3-cr indication 
for a color change of background galaxies found around two 
nearby spiral galaxies. Probing dust reddening induced by 
galaxy halos has not been revisited since then, despite the 
dramatic improvement in data quality and sample size. 

In this work we investigate simultaneously gravita- 
tional lensing and dust reddening of background quasars 
by foreground galaxies and associated large scale structure. 
Our observational results primarily make use of the angular 
cross-correlation between the brightness of quasars and 
the projected density of galaxies. We first recover and 
improve upon previous measurements of the magnification 
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of quasar images by gravitational lensing due to galaxies 
located nearby the light paths using a sample of 85,000 
quasars behind 20 million galaxies d erived from the Sloan 
Di gital Sky Survey (fVork et aljlioool : SDSS) Data Release 
3 ( Abazaiian et al] 2005h . Second, this large sample - 
together with high accuracy photometry in five optical pass 
bands - allows us to detect the presence of dust in the in- 
tervening space and explore its distribution and properties. 
This allows us to study the properties of intergalactic dust 
and provides a way of inferring the abundance of dust in 
the Universe. 

We introduce the formalism of brightness-density cor- 
relations in §2. Data processing and measurements are pre- 
sented in §3 & §4. The astrophysical results are given in §5 
and we summarize in §6. When needed we use Qm ~ 0.3, 
f^A = 0.7 and Ho = 100 h kms~^Mpc"\ 



2 FORMALISM OF BRIGHTNESS-DENSITY 
CORRELATIONS 

Let us begin by considering a background source angularly 
close to a foreground galaxy. This galaxy acts as a gravi- 
tational lens, magnifying the source flux and giving rise to 
extinction and reddening due to the presence of dust sur- 
rounding it. The observed flux is then modifled by the com- 
bination of these effects according to 



(5) 



fol 



fo Me 



(1) 



where fi is the gravitational magnification and tx is the op- 
tical depth for dust extinction at an observed wavelength A. 
The corresponding magnitude shift is 

2.5 

5mA = -2.51og^ + ^ ta . (2) 



In 10 



When fi departs weakly from unity, we can re-express this 
relation as 



Smx ^ 1.08 (ta — Sjj.) 



(3) 



where (5/i = 1 — /i. Thus, magnification and extinction com- 
pete in changing the brightness of background sources. Dust 
extinction is in general wavelength dependent while magni- 
fication is achromatic, so the two effects can, in principle, be 
separated using multi-color data. Below we show how corre- 
lations between the density of foreground galaxies and the 
magnitude of background sources allow us to constrain the 
galaxy-mas and galaxy-dust correlation functions. 



2.1 The galaxy-mass correlation probed with 
magnification 

The galaxy-mass correlation is a powerful tool to probe the 
connection between matter and galaxies, providing direct 
constraints on the dark matter distribution and models of 
galaxy formation. To introduce it, we first define the galaxy 
and mass overdensities: 



5g(x) = -^-l and 5™(x) = 



(4) 



where ng and p are the density of galaxies and matter. The 
galaxy-mass correlation is then defined by 



This cross-correlation can be related to the projected surface 
density of galaxies: 

(5g(<^)E(0 + 0)) = (E(e)) 

= P I egm(\/e^T^) dx (6) 

The first relation indicates that the galaxy-mass correlation 
is equal to the mean mass profile around galaxies, at a given 
separation Q. The second relation is simply a projection of 
the 3-dimentional galaxy-mass correlation introduced above 
in Equation [5] and where — (P' + . 

The mass surface density E can be probed with grav- 
itational lensing. In the weak lensing regime, it is straight- 
forwardly related to the observable magnification, according 
to 



(5/i ~ k/2 = E/Ec 



(7) 



Here k is the lensing convergence and the critical mass sur 
face density is given by 

47r G Da Das 
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(8) 



where Di,Ds and Dis axe respectively the angular diameter 
distances to the lens, the source and between the lens and 
the source. 

As indicated in Equation [3] magnification will affect 
the brightness of background sources and induce a corre- 
lation between the density of foreground galaxies and the 
magnitude of background sources. In order to understand 
the impact on observable quantities, let us consider a given 
area of the sky and let No(m) be the intrinsic magnitude 
distribution of some sources. The photons originating from 
these sources may be defiected by gravitational lensing and 
magnification affects their magnitude distribution such that 



N(m) oc No(m ~ Smi„a) 



(9) 



where Srriina ~ —2.5 log fj, is the induced magnitude shift. 
This leads to an observable mean magnitude shift; 



5mobs = (m) - (mo) 



(10) 



It should be noted that for a magnitude-limited sample of 
sources, the mean magnitude shift induced by a population 
of foreground galaxies, SrUind, differs from the observable 
mean magnitude shift of the individual sources, SrUobs ■ The 
difference between the two depends on the shape of the 
source magnitude distribution N(m) and the limiting mag- 
nitude rrii of the sample. In the case where the induced 
magnitude shift 5m is small compared to the limiting mag- 
nitude of the sample, the difference between the observed 
and induced magnitude shift can be linearized in Srriind and 
we have 



Srrio 



Cs X 5mi„d , 



(11) 



where the coefficient Cs depends on the shape of the mag- 
nitude distribution and the limiting magnitude. In the con- 
sidered limit, it is given by 



(12) 



(see the derivation in the appendix) . If the sample of sources 
is not magnitude-limited, dN/dm(m^) — and Cs = 1, i.e. 
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the measured magnitude shift equals the induced one. For 
magnitude-limited samples, in general we have Cs < 1, i.e. 
the observable magnitude shift is smaller than the induced 
one. It is important to note that for power-law luminosity 
distributions Cs — 0. Therefore, while the brightness of each 
object changes by 5m, the mean magnitude of the sample 
remains unchanged. This is due to the inclusion of sources 
which become brighter than the limiting magnitude. For 
SDSS quasars in the g-hand, down to a limiting magnitude 
of = 21, we find Cs ~ 0.25. 

Let us suppose that a population of background quasars 
and foreground galaxies are sufficiently well separated in 
redshift so that physical correlations between them can be 
neglected. Gravitational lensing will give rise to an apparent 
correlation between quasar brightness and the proximity to 
a galaxy overdensity, which we can write as 



{Smobs)ie) = {Smob4(t>) Sg{(t> + 0)) 

~ Cs {5mind{(t)) 5i,((t> + 6)) 



(13) 



where Sniots = m — (m) is the quasar magnitude fluctuation 
and (5g is the foreground galaxy overdensity at a given an- 
gular distance 9 from a background source. If magnification 
effects are in the weak regime, i.e. E/Ecrit «C 1, the above 
relation reads 



. -1.08 CsP^. 



(14) 



On scales below a few arcminutes, a non-linear treatment 
of the magnification as a function of the de nsity contrast 
is nee ded. Such calculations are presented in lMenard et al.l 
where the authors show that non-linear effects can 
increase the amplitude of the magnification by about 30% 
at a scale of one arcminute. 

For populations of quasars and galaxies selected over 
an appreciable redshift range, the above correlation can be 
expressed in terms of the galaxy-dark matter cross-power 
spectrum using the Lim bers' approximation . Following the 
formalism introduced bv iBartelmannI (|l995h and the nota- 
tion laid out in I Jain. Scranton fc ShethI l|2003 ). we can write 



dX / fcdfc/C(fc,e,x)Pgm(fc,x) 



(15) 



where is the cosmological matter density relative 

to critical, x is the comoving distance, /C is the lensing 
kernel, and Pgm{k,x) is the galaxy-dark matter cross-power 
spectrum. The quantity (5mobs)g(^) therefore probes the 
magnification due to individual galaxy halos on small scales 
and on large scales it constrains the large scale distribution 
of matter in the Universe. 



where 5d(x) is the dust density fluctuation. This cross- 
correlation can be related to the projected dust surface den- 
sity of galaxies: 



{5,{cp) Ed(0 + e)) 



(Sd(e)) 



(17) 



At optical wavelengths, dust extinguishes and reddens 
the light of background sources. The galaxy-dust correlation 
can then be probed by measuring the cross-correlation be- 
tween the colors of background sources and the distribution 
of foreground matter. 

Let us define the magnitude shift of a background source 
measured at a wavelength Aq with respect to the mean mag- 
nitude of the sample. 



5ma{<t)) = ma{<t)) - {rria) 



(18) 



Expanding the formalism introduced in the previous section 
to dust reddening effects, and defining the color excess 



Ea0 = E{Xa — Xfi) = 



Siria — Snip , 

1.08 [r(A„) - r(A;3)] 



(19) 



we introduce the quasar color-galaxy density correlation, 
which can be used to probe the galaxy-dust correlation: 



(E<,^(0)5g(0 + e); 



1.' 



(20) 



If the foreground galaxy and background source populations 
are well separated in redshift, the above correlation provides 
us with information on the mean reddening and therefore 
the mean amount dust around galaxies. The information at 
different wavelengths can be used to probe the shape of the 
extinction curve. For scales on the order of a typical galaxy 
virial radius and smaller, the signal is expected to be domi- 
nated by single galaxies and provides with the average dust 
density profile around galaxies. On larger scales it provides 
information as to the large scale distribution of dust in the 
Universe. 

For a given extinction curve the above quantity can be 
used to infer the mean dust extinction profile around galax- 
ies {A\){8), or similarly the mean optical depth for dust 
extinction {T\)g{9). 

It should be noted that not only extinction but also 
reddening effects can cause quasars to drop out of the 
selection criteria as their identification is a function of 
colors. In the latter case, the measured reddening excess 
turns out to be lower than the true value and the relation 
between observed and induced color change needs to be 
quantified. In the present analysis, reddening effects are suf- 
ficiently small for this to be negligible (see discussion below) . 

The measurement of these galaxy-mass (Equation [S| 
and galaxy-dust (Equation I20|) correlations will be the focus 
of the remainder of this paper. 



2.2 Dust extinction and reddening 

Statistical properties of the distribution of dust around 
galaxies can be constrained by the galaxy-dust correlation 
function: 

C<,d(r) = (<5g(a.)5d(x + r)) . (16) 



3 ANALYSIS 
3.1 The data 

The data set consists of galaxy and quasar catalogues and is 
drawn from the third SDSS data release (DR3; ?). The sur- 
vey provides images in five broad optical bands {u,g,r,i,z 
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Figure 1. Normalized rcdsiiift distributions of tlie galaxies (solid 
blue line) and photometric QSOs (solid red line). The dashed 
line shows the galaxy redshift distribution weighted by the corre- 
sponding lensing efficiency. 



iFukugita et~aLl 1 19961 : [Smith et al.] l2002h . Before masking, 
this set covers roughly 5300 square degrees, the majority of 
which is located around the North Galactic Cap. To reduce 
systematic errors in the photometric data, we impose a see- 
ing li mit of 1"4 and a Galac tic extinction limit of 0.2 in the r 
band l|Scranton et al.|[2002l ). We also apply a mask blocking 
a one-arcminute radius around bright galaxies (r < 16) and 
stars with saturated centers to prevent changes in the source 
number densities due to local fluctuations induced by e rrors 
in sky brightness subtraction l|Mandelbaum et"^ ] [2005l ). Al- 
together, the masks reduce our total area to ~ 3800 square 
degrees. 

With these cuts, we can reliably perform star/galaxy 
separation using B ayesian methods to i — 21 
jScranton et al.l 120021 ') . yielding a galaxy sample whose 
density is independent of local variations in seeing. Galactic 
extinction, sky brightness or stellar density after masking. 
This yields about 24 million galaxies between 17 < i < 21 
at a density of approxim ately 1.8 galaxies per square 
arcminute. The Scranton c t al.l (|2005l l paper used a sample 
that was instead r-band limited, which accounts for the 
difference in the galaxy sample sizes. 

For galaxy magnitudes, we use composite model mag- 
nitudes, constructed from the expMag, devMag and fracDev 
parameters in the SDSS database. This provides a more ro- 
bust estimate of galaxy flux at faint magnitudes than the 
Petrosian magnitudes used in the SDSS spectroscopic sam- 
ple, which can have strong variations due to local seeing for 
r > 18. For quasars, we consider the flux within the psf pro- 
file, designated as psf Mag. All magnitudes are de-reddened 
to correct for Galactic extinction before applying the various 
magnitude cuts. 

The quasar data set was generated using the kernel den- 



sity es timation (KDE) method described in [Richards et ahl 
(| 2004| ) applied to the DR3 data set. The KDE method is 
an extension of the traditional color selection technique for 
identifying quasars. Two training sets, one for stars and 
one for quasars are prepared and the colors for each ob- 
ject are compared to those of the two training sets using 
a 4D Euclidean distance. The objects are then classified 
as either quasar or star according to a larger probability 
of membership. This technique allows a clean separation of 
relatively low redshift {z 5C 2.5) quasars from the stellar 
locus, producing a catalog of 225,000 quasars down to a 
limiting magnitude of g = 21 with efficiency and complete- 
ness grea ter than the SDSS spectroscopic target selc tion al- 
gorithm (iRichards et al.l |2002| : [Blanton et al]|2003l ). After 
masking, the total quasar population is reduced to 195,000. 

El 

In addition to finding quasars, we applied photometric 
redshift techniques ( Wcinstcin et al. 2004) to exclude low 
redshift quasars which might be physically associated with 
our foreground sample. Quasar photometric redshifts are 
driven by the broad spectral emission features in their spec- 
tra, resulting in photometric redshift likelihoods that can 
have multiple peaks as well as strong asymmetries around 
the most likely redshift. Rather than estimating a Gaussian 
redshift error, we use an upper and lower redshift bounds 
for a specified likelihood. To minimize the overlap with the 
galaxies in the redshift space, we require the upper and 
lower bounds to be within the range 1 < z < 2.2. The 
corresponding distribution of quasar redshift probabilities is 
shown in Figure [T] with the red line. For simplicity we treat 
the redshift p.d.f. of each quasar as a flat distribution be- 
tween the upper and lower redshift bounds described above 
and weighted by the likelihood that the redshift was within 
those bounds. This flnal selection criterion reduces the num- 
ber of quasars to 85704. 

We estimate the overall shape of the galaxy red- 
shift distribution bas ed on the CNOC2 luminosity func- 
tions (|Lin et al.lll999l) following the treatment described in 
iDodelson et al.l l|2002l ). This distribution is well fit by the 
expression 



dN 
dz 



^2g-(z/0.187)^ 



(21) 



plotted with the solid blue curve. The mean redshift for the 
galaxy sample is found to be 



{z) ~ 0.36 



(22) 



For lensing purposes it is useful to compute the galaxy red- 
shift distribution weighted by the expected lensing efficiency. 
We have computed this quantity using the quasar redshift 
distribution and we show it using the blue dashed line in 
the figure. The effective redshift weighted by the lensing ef- 
ficiency is given by 



Jdz,dzQZ,i^£Lj^-l^,^,,Q) 
Jdz.dzQ M.j:-J^{z„zq) 



(23) 



It is found to be 2; ~ 0.38. 



^ Upcoming analyses will m ake use of the larger a nd deeper SDSS 
photometric quasar catalog ijRichards et al.ll2007^ . increasing the 
number of objects by an appreciable amount. 
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Figure 2. Left: number counts of our selected sample of photometric quasars with 2 > 1 as a function of magnitude. The solid line 
shows the counts in the g-band, where the quasar sample is magnitude limited. Right: color distributions of the quasars, showing that 
the limiting quasar colors are significantly greater than the modes of the color distributions. Note that for better display purposes, we 
show only objects with colors lower than 0.8 



Figure [2] shows the quasar magnitude and color dis- 
tributions. As can be seen, the quasar sample is magnitude 
Umited in g but not color limited. While quasars are selected 
both in magnitude and color space, the modes of the color 
distributions are far from the limiting colors. As a result, a 
measured color change can be directly used as an estimate 
of the intrinsic color change of the population. Such a prop- 
erty is not satisfied for magnitude changes (see section 12.21 
for more details). 



3.2 Measurement 

We measure the density of galaxies (taking into account 
missing area due cuts on the local seeing variations, bright 
stars, etc.) and compute its correlation, Wa, with the mag- 
nitude of background quasars in the band a, as a function 
of angular separation: 

wc,{e) = {5mc,(4>) 5s{ct, + 6)) . (24) 

On small scales (< 0.01°), the signal contributions from 
magnification and dust extinction are expected to be domi- 
nated by Poisson noise. For larger scales, measurements from 
different angular bins become significantly correlated, as one 
would expect since different quasars will be sampling the 
same local population of galaxies. Since we are looking for 
very small variations in the quasar magnitude and galaxy 
density, the photometric calibration across the survey needs 
to be highly homogeneous. To avoid possible drift in the 
SDSS calibration over time, we measure the QSO-brightness 
galaxy-overdensity correlation separately on each stripe of 



the survey, and then take the average over all stripe-based 
estimators 0: 

^.W = T7^ E "'^'^W' (25) 

JVstripo — ; 

I— 1 

where Wa^ is the QSO-magnitude galaxy-overdensity 
correlation measured for stripe i using the filter a. The 
associated cost is to reduce the sensitivity of the estimator 
to power arising from one direction, i.e., the scanning direc- 
tion, which is equivalent to reducing the size of the sample 
when measuring the large-scale power. Our strongest signal 
is on angular scales smaller than the width of a single stripe 
(2.5°) so the extent of this effect should be minimal. 

In order to measure reddening induced correlations, we 
compute the difference between two stripe-based estimators 
at different wavelengths, 

w^p{e) = wc.{e)-wp{e). (26) 

To estimate the error on the above estimators we use 100 
bootstrap samples of the quasar catalog with which we re- 
measure the above quantities and estimate their dispersion. 

^ Adjacent stripes of SDSS data may be taken on nights with sig- 
nificant observational lag. As such, the photometric zero points 
across stripes are not always calibrated to the precision required 
for this measurement. In principle, this implies that we should 
only use data from single scans over one night. However, in prac- 
tice, we find that stripe-by-stripe treatment yields a sufficent pho- 
tometric zero point accuracy to avoid systematic effects on the 
angular scale of the SDS S stripe. In the future, usin g the photo- 
metric iiber-calibration ijPadmanabhan et ahlEoOSh might solve 
some of the above issues. 
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Figure 3. The anti-correlation between quasar magnitude and foreground galaxy overdensity due to magnification, as a function of scale. 
The light blue and orange data points show the measurements in the g and i bands and their systematic offset suggests the existence of 
dust extinction. We illustrate its effect in the inset where we show, for the five SDSS passbands, the wavelength dependence of the signal 
on small scales. For reference the physical scale at the mean redshift of the galaxy sample is shown in the top axis. 



4 RESULTS 

Following the above procedure we measure the correlation 
between the observed QSO brightness excess and galaxy 
density WaiO) ('Eauationl24[l. where a denotes one of the five 
SDSS pass bands, with the errors estimated from bootstrap 
resampling. The results shown in Figure |3] are examples for 
the g and i bands. We observe that quasar magnitude shifts 
and galaxy overdensities are anticorrelated. In other words, 
quasars appear to be brighter when closer to galaxies, which 
implies a dominance of magnification over extinction effects 
for the main sample of SDSS galaxies with i < 21. At the 
same time, the systematically larger amplitude in the red- 
der band indicates the presence of a wavelength-dependent 
effect which could be attributed to dust extinction. 

Before interpreting these signals, we note that both 
magnification and reddening effects are observed over a wide 
range of scales, from about 5"to a few degrees. We have 
therefore expanded the angular ran ge used in the detection 
of magnification by IScranton et al.l ((2005, ) to both smaller 



and larger separations. Using the mean galaxy redshift com- 
puted in Equation [22] an angular scale of one arcminute cor- 
responds to roughly 150 kpc. The above angular range 
thus corresponds to physical scales from ~ 20 /i~^kpc to 
~ 30 /i~^Mpc, as indicated in the top abscissa of the figure. 

In the inset of Figure [31 we show the wavelength de- 
pendence of the signal on small-scale with measurements 
of Wa{d) with the u,g,r,i and z filters. We can observe a 
continuous trend as a function of wavelength: the signal is 
systematically stronger in redder bands, consistent with red- 
dening by dust. This also shows that probing magnification 
requires a correction for the effects of dust extinction. 

Systematics Tests 

Based on the same quasar sample and a slightly b r ighter 
galaxy sample (selected with r < 21) IScranton et al.l (|2005l ) 
used a density-based estimator and showed that the domi- 
nant magnification signal of the quasar-galaxy correlation 
follows the expected dependence as a function of quasar 
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Figure 4. Correlation between QSO reddening and galaxy overdensity as a function of angular scale. Note that the four independent 
colors are taken from adjacent passbands and do not maximize the signal-to-noise ratio (see Figure|6]for such a quantity). 



magnitude. No such behavior could be detected when us- 
ing stars instead of quasars as a control sample. 

The amplitude of dust reddening depends only on the 
properties and amount of dust around the foreground galax- 
ies selected for the cross-correlation. We can recover a similar 
wavelength-dependent extinction as in Figure |3] or the mean 
reddening signal presented in Figure |4] by using subsamples 
of quasars in different magnitude ranges. As a further san- 
ity check, we also replaced our quasars with stars selected to 
have the same magnitude and spatial distribution on the sky. 
Measuring the cross-correlation between star brightness and 
galaxy overdensity with this sample, we find no appreciable 
color excess for galaxies with i < 20.5 and stars selected 
in various magnitude rangeslf] Finally, we have investigated 
the dependence of the signal as a function of Galactic red- 
dening: by splitting the dataset into tw o regim es of Galactic 
extinction given by the 'Schleg el et alj (|l998h map, we did 
not detect any significant change in our signal. 



WafiiO), where a and (3 indicate two different pass bands. 
We estimate the errors by computing the color covariance 
matrix from bootstrap resampling. Note that the errors on 
colors are smaller than the errors on brightness changes. 
As Figure |4] shows, quasar colors and galaxy overdensities 
are positively correlated, i.e. quasars appear to be redder 
when closer to high concentrations of foreground galaxies. 
The reddening effects are detected for all color combina- 
tions, from 6 ~ 0.1" to about 2°, corresponding to phys- 
ical scales ranging from 50 /i~'^kpc to about 10 /i~'^Mpc. 
The measurement probes galactic radii well beyond the size 
of galactic disks. The amplitude of the effect is stronger in 
bluer bands. We find that a background source whose light 
passes at around 20 /i~^kpc from a foreground galaxy (se- 
lected with i < 21) will be, on average, redder by E{g — i) ~ 
E{B - V) ~ 0.01 mag. 



4.1 Reddening 

4.1.1 Scale dependence 

We isolate the reddening effects by measuring the correla- 
tion between QSO color and foreground galaxy overdensity. 



^ We have detected, however, some excess reddening for stars 
when fainter galaxies (20.5 < i < 21) are used. This effect is 
likely to be attributed to a small contamination of faint galax- 
ies in the stellar sample at faint magnitudes where star-galaxy 
classification is incomplete. This is expected to occur at i ~ 21. 
The presence of galaxies in the star sample may give rise to a 
similar correlation between source reddening and galaxy over- 
density, since galaxy clustering gives rise to an excess of redder 
galaxies in overdense regions, which produces a similar signal that 
mimics the reddening by dust. We have measured the reddening- 
clustering correlation of SDSS galaxies with i > 20.5 and found 
that a contamination of galaxies at about 10% could explain the 
reddening signal seen around stars. 



4.1.2 Wavelength dependence 

The five SDSS filters allow us to constrain the shape of 
A(A), the extinction curve of the dust associated with the 
galaxies, through four independent colors. We measure the 
Wap correlations for two angular bins: O.U < 9 < 0.80 ar- 
cmin and 0.8 < 6 < 8.0 arcmin, which correspond to ef- 
fective projected radii of about 20 < rp < 100 /i~^kpc and 
100/i"^kpc <rp < 1/i^^Mpc. In Figure E] we show the cor- 
responding color excesses with respect to the r band. We 
compare these reddening measurements to the standard ex- 
tinction curve by fitting these data points wit h the func- 
tional form of the extinction curve provided bv lO'Donnelll 
(jl994l ). Such extinction curves are usually characterized by 
the parameter Rv = Ay / E{B — V), which characterizes 
the slope of the extinction curve. The coefficient Ay quan- 
tifies the amount of dust through its extinction in the V 
band. The best fit for Ay and Rv is shown with the blue 
curve. On small scales our measured reddening corresponds 
to Av = (1.3 ± 0.1) X 10"^ mag and Rv = 11.3 ± 7.5, i.e. 
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(28) 



This quantity is proportional to the dust surface density 
and therefore provides us with direct constraints on the 
spherically average distribution of dust around galaxies. 
The angular dependence is similar to average mass profiles 
arou nd galaxies constrai n ed from galaxy-ga l axy l ensing 

([20051)) and 
(2005). This 



(e.g. ISheldon et"al] (|2004h : iMandelbaum et al 



magnification measurements bv lScranton et al 



result has a number of implications regarding the amount 
and nature of the dust in galaxy halos. It will be discussed 
below in section [5] 

In order to convert reddening into extinction, we choose 
the value Rv ~ 3.1, corresponding to standard interstellar 
dust in our Galactic disc and in agreement with the con- 
straints obtained above. In Figure [6] we plot the mean ex- 
tinction profile around the galaxies in our sample. It can be 
written as 



(Av)(6l) = (2.4 ±0.2) X 10" 



(29) 



Figure 5. Mean color excess with respect to the r-band measured 
in two angular bins. The blue curves sh ow the best fit extinction 
curve parametrized bv lO'Donneiil l ll994) . The dashed ereen curves 
are fits with Ry = 3.1. Note that Ry is poorly constrained on 
small scales. 



the slope of the extinction curve is not well constrained. On 
large scales however we obtain a better accuracy: 



(Av)(r-p) = (4.14±0.19)xl0" 



Ry = 3.9 ±2.6. 



(27) 



and Ay = (1.2 ± 0.3) x 10"^ mag. The green dashed curve 
shows the best fit for Ay when Ry is assumed to be 3.1, the 
standard value for dust in the disk of our Galaxy. Within the 
errors our results are consistent with standard interstellar 
dust. 



4.2 Dust extinction 

Given our measurements of the reddening due to foreground 
galaxies, we can now make estimates of the average dust 
extinction profile of the galaxies in our sample. In addition 
to constraining the amount and spatial distribution of dust 
on large scales around galaxies, this quantity is of interest 
to quantify the intrinsic brightness of background sources, 
which is the main goal of Type la supernovae measurements 
aimed at constraining dark energy. 

As indicated in ^2.21 the correlation between quasar 
color and galaxy overdensity is an estimator of the mean 
reddening induced by galaxies. In order to maximize the 
signal-to-noise ratio of the reddening detection, we measure 
the quantity Wap{d) for which a and 13 are taken to be the g 
and i pass bands. We avoid using the u and z bands for which 
the photometric errors are substantially larger. In addition, 
the latter sufi'ers from photometric contamination due to sky 
emission lines. We show the mean E{g — i) color excess of 
quasars as a function of scale in Figure |6] A best fit power- 
law gives 



lOOft-ikpc 



(30) 



where Ay is the observed V-hand extinction. For an extinc- 
tion curve characterized by Rv = 3.1, we have Aa oc X~^'^ 
in the visible range. We point out that the amount of red- 
dening profile shown in Figure [6] is significantly lower than 
contributions expected from dwarf galaxies. For example, 
the average dust reddening induced by the LMC, located 
at ab out 50 kp c from our Galaxy, is E(B — V) ~ 0.075 mag 
(Schl egel et al.i ri998V Such a value is about an order of mag- 
nitude larger than eq. [28] and shows that satellite galaxies 
are expected to dominate the amount of dust reddening for 
individual lines-of-sight intercepting galaxy halos. 



4.3 Magnification 

The correlation between quasar magnitude shift and 
galaxy overdensity allows us to constrain a combination of 
magnification and dust extinction. Having estimated the 
wavelength-dependent extinction from reddening measure- 
ments in the previous section, we can now use this informa- 
tion and isolate the magnification effects. Our goal here is 
to estimate the signal well enough to allow for a compari- 
son with similar measurements from galaxy-galaxy lensing 
available in the literature. 

For simplicity we start by describing the angular cor- 
relation between quasar brightness and galaxy overdensity 
by a simple power-law with index —0.8. As above, we use a 
value of Ry — 3.1 (or Aa cx A^^'^ in the visible range) to 
convert reddening into extinction. The mean quasar bright- 
ness change as a function of angular separation from galaxies 
then reads 



(5mob.)(6') 



Sm/j, ± Sm-T- 



1' 



(31) 



where (5m^ is the achromatic magnitude change due to mag- 
nification and SttIty is the change due to dust extinction in 
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Figure 6. Mean observer-frame E(g-i) reddening profile around galaxies with i < 21, as a function of impact parameter in angular scales 
(lower axis) and effective projected radius (upper axis). Note that we have E(g — j) ~ E(i3 — V). The right axis shows the corresponding 
observer-frame extinction in the V band using the standard interstellar value Ry = 3.1. 



the V-band. As shown above, we have &m,Ty — +2.4 x 10 
(Eq.[29ll, which gives 5m^ ~ -6.7 x 10"^ 

At sufficiently large wavelength, dust extinction be- 
comes negligible and the brightness change is only due to 
magnification. In this limit, we find 



(32) 



~ 2.5 X 10-2 



where we have used the conversion between observed and 
induced magnitude shift (Equation I14|) and Cs ~ 0.25 for 
SDSS quasars with g <21 (Equation I12p . This implies that 
the mean magnification excess found at one arcminute (or 
~100 kpc) is abou t 2%. T his result is in agreement with 
that of lScranton et af] (|2005l ) where the authors inferred the 
magnification from the observed correlation between quasar 
and galaxy densities on the sky, i.e. a different estimator. 
Our analysis also allows us to measure magnification effects 
on smaller scales. On the smallest scale we can probe, i.e. 
5" (or ~15 h.-^ kpc), we find S^i ~ 15%. 

In the weak lensing regime, the magnification is related 
to the mean mass surface density by 



1 + S^i 



(33) 



The inferred mean mass surface density of our galaxy sample 
is shown with the blue data points in Figure [T] A best fit 
power-law distribution gives 



1' 



(34) 



with A = 30.6 ± 3.4 /iMq pc or using our effective pro- 
jected scale: 

-0.8 



E(rp) ~ A' 



Ift-iMpc 



(35) 



with A' — 8.1 ± 0.9 ft Af0 pc~'^. These results provide an 
estimate of galaxy density profiles from gravitational mag- 
nification corrected for dust extinction effects. While a lin- 
ear relation between magnification and density contrast is a 
good approximation on large scales, higher-order corrections 
become signi ficant on scales smalle r than a few arcminutes. 
As shown by iMenard et al] (|2003l ). using Eq. 1331 results in 
overestimating the projected mass E by about 15 to 25% on 
scales ranging from 1 to 0.1 arcminute. The first few points 
shown in Figure [7] have not been corrected for this effect. 
As shown below, on larger scales the mass estimate from 
magnification is in good agreement with shear-based mass 
estimates. 



4-3.1 Comparison with shear measurements 

So far the galaxy-mass cross-correlation has mostly 
been accessible through galaxy-galaxy lensing measure- 
ments which estimate the mean tangential shear of 
background galaxies for a given s ample of fore g round 
lenses (e.g. iHoekstra et all I2OO2I , ISheldon et all |2004| , 
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Figure 7. The mean surface density of galaxies (with i < 21) measured through the magnification of background quasars and corrected 
for dust extinction (blue points). In comparison we s how the mean s urface density of a sample of ^ L* galaxies at 2 ~ 0.1 obtained from 
the gravitational shear of background galaxies from IS heldon et al. I ||2004|) . Non-linear magnification effects have not been included and 
result in an overestimation of the mass on the smallest scales. 



iMandelbaum et all [20051 . iLeauthaud et all 12007! ). The tan- 
gential shear (74) azimuthally averaged over a thin annulus 
at projected radius R from a lens galaxy is directly related 
to the projected surface mass density of the lens within the 
aperture, 



7t 



AE(il) 



where 

AE(i?) = E(< R) - T,{R) , 



(36) 



(37) 



E(< R) is the mean surface density within radius R, and 
T,( R) is the azimuthally a veraged surface dens ity at radius 
R (jMiralda-Escudelll99ll : iFahlman et al.lll994l ). Therefore, 
shear measurements constrain AE whereas magnification is 
a direct estimate of E. 

In order to compare magnification and shear measure- 
ments, it is interesting to note that, in the case of power- law 
mass profiles, with E(rp) oc r~", we have 



AE(rp) 



■ S(r,) , 



(38) 



The two observables E and AE are therefore equivalent 
for isothermal profiles. For an angular dependence follow- 
ing r-°-^, we have AE = E x / with / ~ 0.66. This 
value goes down to 0.54 for an index of —0.7. Using galaxy- 
galaxy measurements based on the SDSS, ISheldon et al.l 
l|2004h found AE(i?) = A' {vp/l h'^Mpcy"' with A' = 
(3.8 ± OA)hMQpc'^ and a' = 0.76 ± 0.05 for a sample 
of spectroscopically identified lenses with (z) ~ 0.1 and 



(L) ~ L* . Considering for simplicity a slope of —0.8 (con- 
sistent with their constraint), their results translate into 
E(i?) = (5.7 ± 0.6) (rp/l/i-iMpc)"°-%A'/Qpc-2. Their 
scaled-measurements (using Equation 138 |l are shown in Fig- 
ure[7l Under these assumptions, the magnification and shear 
estimators appear to be in very good agreement. A more ac- 
curate comparison between the two would require account- 
ing for the differing luminosity distributions between the 
two lensing samples as well possible redshift evolution (our 
lenses comprising most of the sources for the shear-based 
estimators). However, it illustrates that both methods offer 
comparable degree of measurement precision and dynamic 
range while affected by significantly different systematic ef- 
fects. 



5 IMPLICATIONS 

Our detection of the change in apparent magnitude of 
distant quasars has allowed us to quantify the magnification 
and reddening of background sources as their light rays 
pass in the vicinity of foreground galaxies. While the am- 
plitude of gravi tational lensing effe c ts has been known both 
from theory (iBartelmannI 1 19951: iMenard fc BartelmannI 
l2002l: IJain Scrantoii fc ShethI l2003h and observations 
I Scranton et al. 2005h . the expected amplitude of dust 
reddening and extinction effects was largely unconstrained. 
Our study has shown that, in the visible bands, dust extinc- 
tion occurs at a level comparable to that of the observed 
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brightening due to magnification. Hence, studies aimed at 
predicting observable magnitude changes of background 
sources (quasars, galaxies, supernovae, etc.) which include 
only gravitational lensing effects are incomplete. In the 
visible range, dust extinction effects cannot be neglected 
and must be included alongside with magnification to 
properly account for the effects of passing through large 
scale structure. 

Below, in order to simplify the discussion we will 
assume the dust in galactic halos to be described by SMC 
type dust. This choice is motivated by a number of results: 
(i) certain low-ionisation absorbers such as Mgll are known 
to inhabit the halo of ~ L* galaxies. They are found on 
scales reaching up to about 100 kpc around galaxies 
l|Zibetti et al.|[2007l ). Recently several a uthors studied their 
extinction properties and dust content (jMenard et al.ll2005l : 
[Kharo ot al. 2005; York ot al. 200i). They found that their 
average extinction curve is similar to that of the SMC, i.e. 
does not show the 0.2 ^.m bump seen in the Milky way 
extinction curve (for systems with z > 0.9 where the feature 
enters the visible window), (ii) In addition, it is known 
that only a small fraction of high redshift galaxies show an 
extinction curve with the 0.2ii.m bump. As discussed below, 
considering a Milky- Way type dust changes our results only 
by a factor two. 

In order to characterize the population of galaxies re- 
sponsible for most of the intergalactic dust, we first assume 
that the amount of dust in halos is, on average, proportional 
to the metallicity and luminosity of galaxies. Under this as- 
sumption, the dominant contribution of the reddening signal 
is expected to originate from galaxies with an effective lu- 
minosity: 



(^(L)Z(L)LdL 



'i: 



dL0(L)Z(L) (39) 



where <?!>(I') is the galaxy luminosity density given by the 
Schechter function 



<t>{L)dL = <t>* ( — 



-L/L' 



(40) 



with a — —1.1 and Z(L) i s the metallicity-luminosity rela- 
tion l|Tremonti et al.ll2004l '): 



12 



log(g 



-0.185 Ms +5.328 



(41) 



with Mg = —19.5. By integrating Equation |39] between 
L — 10~^ L* and infinity, we find that the most important 
contribution of reddening originates from galaxies with a 
luminosity 



Left ~ 0.45 L" 



(42) 



Given the nurnber density of L* galaxies 
iFukugita fc PeeblesI |2006| ) and the above luminosity 
function, galaxies with L = Loft are expected to have a 
comoving number density of n ~ 0.037 Mpc~^. 

Below we discuss a number of implications given by 
the existence of a large-scale distribution of dust around 
galaxies. We note that only certain of the following results 
will make use of the above numbers. 



5.1 Dust mass distributions 

The properties of dust particles giving rise to the observed 
optical extinction and infrared emission in our Galaxy, the 
LMC and the SMC ha ve been modeled by several au - 
thors (e.g. iMathij [ioool and IWeingartner fc Drain3 1200 il l. 
These models provide estimates of the absorption cross sec- 
tion per mass of dust as a function of wavelength _R'ext(A), 
which is used to infer dust mass surface density from 
reddening profiles. For the SMC type dust the model of 
IWeingartner fc Draind |20o3) gives M 



Kext(A\ 



1.54 X lO'^cm^g"^ 



(43) 



Note that Milky Way type dust corresponds to a dust mass 
larger by a factor two at a fixed extinction. 



5.1.1 Evidence for a 
halos 



fuse component of dust in galaxy 



The knowledge of Kcxt allows us to convert the observed red- 
dening into a dust mass surface density. The spatial depen- 
dence of this quantity is shown in Figure[8l We find a column 
density of dust of about 10~^/i M© pc~^ at an impact param- 
eter of 100 kpc. It is interesting to estimate the total 
amount of dust in the halo of the galaxies defined by Equa- 
tion l42l Considering an isothermal sphere mass distribution, 
the virial radius of 0.5 L* galaxies is r„ ~ 110/i~^kpc. The 
mass of dust residing in the halo of such galaxies is given 
by integrating the ratio Av{rp)/ Kcxt{^v) over the area en- 
closed by r„: 

2-K In 10 



Mdust — 



2.5 Kext(Av) 



r 

Jo 



Av(rp) rp drp 



(44) 



Our measurements provide us with the radial dependence 
of Av for r-cff > 20 kpc, i.e. scales greater than galactic 
disks. By defining the dust mass in the halo as M^„!j° = 
Md„st(20/i"^ kpc < reft < r„) and using eg. [29] with SMC 
type dust, we find 



mS ^ 5 X 10' Mq 



(45) 



This estimate shows that a substantial amount of dust exists 
in the halo of ~ L* galaxies. This dust mas s is comparable t o 
that commonly found in galactic disks (Dr aine et al.|[2007l ). 
This immediately shows that the reddening profile shown in 
Figure [S] cannot be explained by a contribution from satel- 
lite galaxies whose luminosities and therefore dust masses is 
negligible compared to a central L* galaxy. Our results sug- 
gest the existence of a diffuse component of dust in galactic 
halos. 

We note that similar arguments can be obtained from 
reddening considerations and do not rely on the specific 
value of Kaxt (Eq. 143 |l . It should also be noted that while 
the total amount of dust appears to be dominated by a dif- 
fuse component, the largest reddening values are still ex- 
pected to originate from lines of sight passing through satel- 
lites. The average reddening value of the LMC, E(B — V) ~ 
0.075 mag, is about one order of magnitude larger than the 
average halo reddening at a radius on order 50 h"^ kpc (see 
Figure [6ll . 

^ We thank Joseph Weingartner for having provided us with this 
value. 
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Figure 8. Correlation between magnification and galaxy over- 
density as a function of scale. On scales smaller than ~ 500 
/i~^kpc, this quantity is a direct estimate of the mean surface 
density of galaxies (with i < 21). 



5.1.2 Dust-to-mass ratio 

Having obtained constraints on both the total mass distri- 
bution from magnification and dust mass from reddening 
(see Figure [S}, we can compare their statistical spatial dis- 
tributions by computing the ratio 



m 



(Sdust(^)) 



(46) 



where Edust and E are the dust mass and total mass surface 
densities. This ratio is plotted as a function of scale in the 
lower panel of Figure |8] and appears to be only weakly scale- 
dependent. We find the dust to mass ratio 



r ~ 1.1 X lO""" . 
for SMC type dust. 



(47) 



5.1.3 The cosmic density of dust 

We now attempt to estimate the cosmic density of dust, 
^dust. To do so we first compute the density of dust 
originating from galaxy halos. Considering ~ L* galax- 
ies with a mass-to-light ratio M„/Lb ~ 250 h Mq/Lq 



jFukugita fc PeeblesI 120041') ■ a light d ensity £b ^ 2 x 
10^ feL^Mpc"^ l|Blanton et all I2OO3I ') and for a dust-to- 
mass ratio r in halos, we can write 



i'dust 



r X (m^/Lb) X Ci 



Pcrit 



2.8 X 10" 



(48) 



Previous attempts to obtain observational constraints on the 
cosmic density of dust have focused on light-weighted es- 
timates, i.e. dust related to disks. For example, using the 
attenuation-inclination r elation for galaxy d iscs and their 
associated central bulges, [Driver 

ilal]|2003) 

quantified the 

mean attenuation of a large sample of galaxies and esti- 
mated fl^uBt — 2 X 10"''. Similarly, theoretical estimates 
have focused on dust associated with cold gas. For example, 
iFukugita & Peebles (2004) estimated the cosmic dust den- 
sity by computing the mean metallicity of galaxies weighted 
by the Schechter luminosity function, used a mass fraction of 
metals into dust grains of Z(dust)/Z = 0.2 multiplied by the 
density parameter in cool gas. They find fidust — 2.5 x 10~^. 
The total cosmic dust density can be estimated by summing 
up the halo and disk contributions. We find that 



f^d 



"dust + "dust 

5.3 x 10"® . 



(49) 



This value is in agreement with the (model-dependent) up- 
per limit obtained by Inoue fc Kam aya (20041 : fidust < 10~^ 
at 2 ~ 0. 3. It is about a factor two larger than the estimate 
given bv lFukugita fc Peebles! l|2004h for galactic disks. 



5.2 The opacity of the Universe 

As our results provide us with an estimate of the mean opti- 
cal depth for dust extinction due to a population of galaxies 
spanning a given redshift range, we can attempt to inte- 
grate T^{z,\) over redshift and estimate the opacity of the 
Universe. 

Assuming that most of the dust in the Universe is asso- 
ciated with galaxies and that most of the optical depth for 
extinction originates from the halo of galaxies with L ~ L^g 
(defined in Equation l42p . the mean dust optical depth in the 
Universe up to a redshift z is then given by 



fix, 



A \ c(l + 2)^ 



1 + 2 



Hiz) 



dz 



(50) 



where a and n are the cross-section and number density 
of galaxies with L ~ I/off- For a dust profile extending up 
to the virial radius, r„, of these galaxies we have ct = 7rr^ 
and fg is the average optical depth within the halo, which 
we define with 20/i~^kpc < r < r„. We evaluate f{\,z) 
(Equation I50|l and present the results in terms of extinc- 
tion in the observer-frame V-band. Integrating Equation 
1501 up to the mean redshift of the galaxy sample used in 
this analysis, we find the value denoted by the blue circle: 
(Av(2 ~ 0.36)) ~ 0.02 mag. We then estimate Av(2) for 
different dust models shown in Figure [9] 

(i) we first use a constant dust density with redshift, rep- 
resented with the dashed line. At high redshift, such an esti- 
mate is expected to represent an upper limit on the allowed 
extinction. 
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Figure 9. The average observer-frame Ay extinction as a 
function of source redshift. The blue circle shows our (model- 
dependent) estimate of the opacity induced by dust in the halo of 
~ L* galaxies at the mean redshift of our sample. The three blue 
curves represent extrapolations to higher redshifts, using different 
dust models: a constant co-moving density (dashed line), a den- 
sity decreasing with the mean metallicity (solid light blue) and a 
suppressed density at high redshift (dark blue). Various observa- 
tional constraints are shown and described in Appendix B. 



(ii) A more realistic estimate should take into account the 
fact that the amount of dust in and around galaxies is red- 
shift dependent as dust is being produced at a rate following 
that of the metals. As mentioned above, useful constraints 
on the evolution of the amount of dust around galaxies come 
from studies of metal absorbers. Recently, I Menard et al.l 
l|2008l ) probed the redshift evolution of the amount of dust 
associated with Mgll absorbers from z — 0.4 to z = 2. They 
found that their dust content follows pdust oc (1 + z)~^'^ , 
which turns out to be similar to the evolution of cosmic star 
density. By taking this redshift dependence into account, we 
obtain a an alternative estimate of the cosmic opacity. The 
corresponding results are shown in Figure |9] with the solid 
light-blue curve. As can be seen, taking this effect into ac- 
count brings a modest change to our previous estimate and 
lowers the total opacity by about a factor two at 2 ~ 2. 

(iii) Finally, to illustrate the range of possible redshift de- 
pendences, we consider a third model where we damp the 
density of dust by an addition factor (1 + z)~^ , shown with 
the dark blue curve. This model is somewhat ad hoc in na- 
ture, but demonstrates that the current measurements and 
limits allow for considerable variation at higher redshift. 

Figure|9]also shows various upper and lower limits (detailled 
in Appendix B) on the opacity as a function of redshift. As 
can be seen, at z ~ 1, Ay values are bounded within about 
a factor 10. 



5.2.1 Implications for supernova experiments 

Type la supernovae provide us with an estimate of luminos- 
ity distances and have been extensively used to constrain 
cosmological parameters, dark energy in particular. Super- 
novae are standardizable candles and become usable as stan- 
dard candles after re-normalizing their brightness for intrin- 
sic brighter-bluer and brighter-slower trends as well as dust 
extinction. We now investigate how the presence of inter- 
galactic dust can affect such constraints. 

The color c of each supernovae is the sum of several 
contributions: c — X^^Ci, where Ci are the intrinsic color, 
dust reddening by the host, dust along the line-of-sight, 
etc. Each of them can be corrected for using an appropri- 
ate "reddening-to-extinction" coefficient The observed 
supernova magnitudes are used as a distance estimator ac- 
cording to 



M + a{s, - 1) -/3ci 



(51) 



where the apparent magnitude mi, the stretch Si and the 
color Ci are derived from the fit to the light curves. The 
parameters a, /3 and the absolute magnitude M are fitted 
by minimizing the residuals in the Hubble diagram. A color 
excess Ci which does not contribute to a significant scat- 
ter will not affect the inferred value of /3. Its reddening-to- 
extinction coefficient will be described by the best-fit /3o and 
might lead to a bias if its intrinsic /3 differs from this value. 
As shown above, the presence of intergalactic dust might 
introduce such an effect. On large scales around galaxies 
we find Rv ~ 4, which corresponds to a value of /3 ~ 5 , 
while SN analyses l ead to /? ~ 2 - 3 (|Conlev et al.l l2007l : 
iKowalski et al.! '2008). This indicates that, at some level, the 
(redshift-dependent) contribution of intergalactic might not 
be properly corrected by using Eq. 1511 The magnitude of 
this effect and its impact on dark energy constraints must 
be investigated and quantified. 



5.3 Dust-to-light ratio 

We can compare the statistical properties of the light (traced 
by galaxies) and dust distributions by computing the ratio 
of the galaxy-reddening cross-correlation to the autocorrela- 
tion of the same galaxies. To do so we compute, as a function 
of scale, the parameter 



(52) 



where we use the g — i color as our estimate of dust redden- 
ing (see section r4.2|) . This quantity is shown in Figure [TOl 
The dust reddening-to- light ratio appears to be weakly scale- 
dependent. On scales greater than about one arcminute, we 
find 



f3(3-^)ie > 1') =i 0.015 mag 



(53) 



The excess seen on the smallest scale might be due to a 
contribution from galactic disks. The dust-to-light ratio is a 
quantity which can be used to constrain models describing 
the transport of dust outside of galaxies. Our results show 
that, on large scales, the dust distribution follows that of 
galaxies. 
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Figure 10. Ratio between tlie reddening-galaxy correlation and 
the galaxy autocorrelation function, indicating that the dust to 
light ratio is weakly scale dependent. 



6 SUMMARY 

Intervening galaxy halos and large-scale structure affect the 
light of background sources through both gravitational lens- 
ing and dust extinction. In this paper we present simulta- 
neous detections of these two effects obtained by measuring 
the cross-correlation between the brightness and colors of 
about 85,000 z > 1 quasars and 20 million z ~ 0.3 galax- 
ies from the SDSS. We find that quasar brightness is cor- 
related with galaxy overdensity, indicating that magnifica- 
tion effects dominate over dust extinction for this sample of 
galaxies. By examining the quasar colors, we observe that 
quasars appear to be redder as the projected galaxy density 
increases. Both effects are detected on scales ranging from 
0.1 arcmin to about 2 degrees, corresponding to projected 
radii of 20 kpc to 20 /i~^Mpc at the mean redshift of the 
galaxy sample. More specifically, we find that 

• galaxies and large-scale structures at z ~ 0.3 in- 
duce an excess reddening to background sources given 
by (E(B-V)> ~ 1.4 X 10"^ (S/l')"" *". On arcminute- 
scales around galaxies, the slope of the extinction curve 
is constrained by Rv = 3.9 ± 2.6, a value consistent 
with that of our Galaxy. Assuming the Galactic value of 
Rv = 3.1, we find a mean extinction profile given by 
(Ay) ~ 2.4 X 10"^ (rp/lOO/i-^kpc)"" *". The ratio of the 
galaxy-reddening to galaxy-galaxy correlation functions is 
found to be roughly scale independent, with an amplitude 
(E(sr - i))/wgg ~ 0.015 mag. 

• The amplitude of dust extinction in the 1^-band is about 
one-third of that of the observed brightening due to magni- 
fication. We estimate the extinction-corrected magnification 
profile and fi nd (ti(e)) ~ .025 (g/l')""'^, consistent with 
the results of lScranton et al.. (2005) and extending those to 
both smaller and larger scales. The average mass surface 
density profile of the galaxies inferred from our measure- 
ments is comparable to galaxy-galaxy lensing estimates. 



• At a separation of 20 /i ^ kpc from a 2: ~ 0.3 galaxy, 
a background source is, on average, magnified by a factor 
fi = 1.15 and reddened by E(B-V)~ 0.01. 

The detection of dust on large-scales has a number of 
implications: 

• The amount of dust found in galactic halos is found 
to be comparable to that in the disk of ~ L* galaxies and 
therefore significantly larger than that of dwarf satellites. 
This implies the existence of a diffuse component of dust in 
halos, predicted by some models of dust halo dynamics but 
heretofore unobserved. 

• Having argued that the dominant contribution of dust 
reddening observed in this analysis is due to galaxies with 
L ~ 0.51/*, we have estimated their halo dust mass to be 
about 5 X 10^ Mq within their virial radius. 

• Including both disk and halo contributions, we find 
Sldust — 5 X 10~^, a value roughly twice that estimated by 
(|Fukugita fc Peebles! |2003 ) for galactic disks. 

• Such an extended distribution of dust around galaxies 
will affect the apparent magnitude of distant sources. We 
have estimated the mean opacity of the Universe due to 
dust in galactic halos. Our model-dependent estimate gives 
Av(z = 0.5) ~ 0.03 mag. Such a value is less constrained 
at higher redshifts due to our limited knowledge of the evo- 
lution of dust with redshift. Considering several models, we 
found Av{z = 1) ~ 0.05 - 0.09 mag. This wiU affect the 
brightness estimates of Type la supernovae at high redshift, 
which require high precision in order to maximize their con- 
straints on cosmological parameters. Dust reddening may 
also induce dispersion in brightness and color. 

The technique presented in this paper provides us with 
a unique probe of the distribution of dust (warm and cold) 
on large scales around galaxies, which is otherwise difficult 
to explore. It opens up the way to studies of the amount 
of circumgalactic dust as a function of galaxy type, lumi- 
nosity and environment which may shed light on the origin 
of the dust. Extending the analysis to UV measurements is 
an important task to obtain better constraints on the dust 
properties. The final remark is that this analysis requires 
only accurate photometric data in several passbands. 
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APPENDIX A: DERIVATION OF Cs 

The observable magnitude shift defined by Equation [10] is 
a function of the shape of the magnitude distribution, the 
limiting magnitude and the induced magnitude shift Sm. In 
the case where the induced magnitude shift 5m is small com- 
pared to the limiting magnitude of the sample, the difference 
between the observed and induced magnitude shift can be 
linearized in Sm. We have 



Amobs = 



(m) - (mo) 

J Ammnim ~ Sm) 



J Ammnijn) 



J dmn{m — Sm) J dmn{m) 



(Al) 



For induced magnitude shifts small compared to unity, we 
can Taylor-expand the above expression to first-order in Sm: 



Amobs — -r- 



5m 



J- j Amn'im) — j Amn'(m)m 



(A2) 



where la ~ J dmn(m) m°' . The above expression can be 
simply written as 



Amobs 
with 

Cs = 1 



5m X Cs . 



1 dN 

dm 



(mi) ^mi - (mo) 



(A3) 



(A4) 



We can verify that in the case of a power-law luminosity 
function, i.e. dN/df oc /" or dN/dm oc a"", the above ex- 
pression gives Cs = which implies that no magnitude shift 
can be observed. 



APPENDIX B: OBSERVATIONAL 
CONSTRAINTS ON THE COSMIC 
TRANSPARENCY 

In this appendix, we detail the values of the cosmic opacity 
upper and lower limits used in Figure|9] Note that the upper 
limit values plotted in the figure reflect the 95% or 99% 
confidence limits for those measurements, while the lower 
limit values are taken from the actual measured values. 



• iMore. Bow, fc Hog"3 l|2008l ) put a virtually assumption- 
free constraint on the opacity of the Universe at low redshift 
using the Tolman test: 



Dl = (1 + 2)' Da, 



(Bl) 



where Dl is the luminosity distance and Da the angular di- 
ameter distance, independent of world model. Any observed 
deviation from this expected relation is taken to be a result 
of extinction along the line of sight. Combining observational 
results from supernovae and baryon acoustic oscillations to 
estimate the change in optical depth from redshift 0.20 to 
0.35, they found that Ar < 0.13 at 95% confidence. As- 
suming no evolution of the dust properties in the redshift 
range < z < 0.25, we can use their result to estimate the 
expected Ay extinction up to 2; ~ 0.35: 



Av^l.08 Arx [l+ir/Io:f] 

where I^ = /J d^ (1 + z)yH{z) 
0.14 mag. 



(B2) 



. This gives Av{z = 0.35) < 



• lAvgoustidis. Verde, fc Jimenej l|2009l ) obtained a con- 
straint on the low-redshift opacity by combining observa- 
tional results from supernovae and the Hubble Key project. 
They simultaneously fitted for H{z) and the luminosity dis- 
tance, where the latter quantity was allowed to be mod- 
ulated extinction. Us ing the Union sample of supernovae 
(|Kowalski et al.ll2008l ) with a mean redshift of z ~ 0.5, they 
found 



Av ~ 1.08 X 2ez 



(B3) 



with e = -0.0ll!5:J]^. This gives Av{z ~ 0.5) < 0.08 mag. at 
99% confidence^ 

• iMortseU fc Goobaij l|2003l ) analyzed the scatter in 
quasar colors as a function of redshift, attributing its ex- 
cess to dust extinction. They reported a 99% upper limit on 
the cosmic opacity : Av(z = 1) < 0.2 mag. 

• iMenard et al.l (I2OO8I ) quantified the mean amount of 
reddening and extinction induced by strong Mgll absorbers, 
i.e. systems usuaUy fou nd within ~ 100 kpc of ~ L* galax- 
ies (jZibetti et al.l [iooj) . They were able to constrain the 
mean reddening {Eb-v) as a function of the rest equivalent 
width Wo of the absorbers and their redshift, in the range 
0.4 < 2 < 2. In addition, as described in section O various 
authors have shown that, on average, the extinction curve 
associated with Mgll absorbers is consistent with that of 
the SMC, i.e. with Rv — 3. The mean extinction induced 
by these systems can then be computed according to 



Av{z) = Rv / AWo / Az — — 
Jo Jo AWoAz 



{Eb-v{Wo,z)). (B4) 



Using the parameterization of the the incidence r ate of Mgll 
absorbers, A^N/AWoAz, given bv iNestor et al.l l|2005h . and 
assuming no evolution in the dust properties in the range 
< 2 < 0.4, we can use these measured values as lower limits 
on the global extinction. We find: Ay (2 = 0.6) > 0.009 mag, 
Av(2 = 1.3) > 0.029 mag and Av(2 = 1.8) > 0.044 mag. 
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